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Abstract 

Topological quantum phase transitions in superconductivity are discussed on two dimensional 
lattices. The main focus is on the Chern number for superconducting states. Each superconduc- 
tivity is characterized by the Chern number, and the quantum phase transition is associated with 
topological changes of the quasiparticle Bloch function in the Brillouin zone. For the supercon- 
ducting case, the Chern number has several equivalent but different topological expressions given 
by vortices, the Dirac monopole and strings. We demonstrate quantum phase transitions by these 
topological quantities both for singlet and triplet cases. 
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The quantum phase transition is a drastic change of a ground state when physical param- 
eters of a system are varied. This is a phase transition at zero temperature with quantum 
fluctuation which raises many interesting physical questions. Recently substantial number of 
works have focused on the topic reflecting its importance. Typical examples include the Mott 
transition in strongly correlated systems and the Hall plateau transition in the quantized 
Hall effect 0, 0, The latter is special in the sense that the ground state is characterized 
by a quantized quantity ( quantized Hall conductance axy )■ As is known today, its quanti- 
zation originates from a topological character of the Hall conductance!^, |^, §]. Then the 
transition is governed by topological objects such as the Chern numbers, vortices and edges 
states[^, 0, ^, 0. Each phase is characterized by the different topological quantum number 
although the symmetry of the state is not necessarily different. It is quite different from 



usual phase transitions and is characteristic to the topological quantum phase transition[10 



There have been several features for this type of the phase transition. Some of them are 
a special form of a selection rule, and the possible stability of the phase which puts spe- 
cial importance on the topological transitions and discriminates from other quantum phase 



transitions H, 10 



Recently, there have been trials to extend the concept of this topological phase transition 
to unconventional singlet superconductivity |n|, |T^, |TB|, |T^. In the discussions, mapping the 



system into the standard quantum Hall system is essential. Here the "spin" Hall conduc- 
tance plays a main role which is given by the Chern numbers. A non trivial example of such 



states is given by time-reversal symmetry breaking superconductivity [0, |T^, [IJ, |T5|, |16 



Leaving its reality apart, the topological phase transition in superconductivity raises impor- 
tant theoretical questions how the topological character restricts transition types. In this 
paper, we take the generalized Bogoliuvov-de Gennes (B-dG) hamiltonian for superconduc- 
tivity to discuss the topological quantum phase transition and demonstrate the topological 
character using several topological expressions. 

Search for spin triplet superconductivity has also long history but several important 
materials have been found recently which raised intense re-investigations of triplet super- 
conductivity. One of the interesting aspects of the triplet state is its rich structure of the 
order parameters which allows the existence of non-trivial topological structures. On this 
point, the topological characters of the chiral p-wave superconductivity is investigated by 



several groupsO, IT^, IT8[. In this paper, we treat this spin triplet case on a two dimen- 
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sional lattice based on the same strategy as the spin singlet case. Then we can characterize 
the triplet state by the Chern number and several topological expressions as well ||19|| . 

Comparing to the usual Quantum Hall effect with many Landau bands, there are only 
two quasiparticle bands of quasiparticles ( and quasiholes ) in the superconductivity. This 
special but simple situation allows us to use the interesting Berry's parameterization. Then 
the gauge fixing to calculate the explicit Chern number also produces an interesting topo- 
logical expression. Using this gauge fixing, we establish a relation between several different 
topological expressions. 

The Chern number To demonstrate the topological transition in superconductivity, we 
investigate systems described by the following B-dG hamiltonian for superconducting quasi- 
particles 

ij icr 

In this letter, we assume the system is translationally invariant as tfj = f^ii — j) and 
Aij = A{i — j). Then the hamiltonian is written as 

' eT(fc) A(fc) 



A*{k) -e^{-k) 



H = ^c\k)h{k)c, h{k) = 

k 

where c^{k) = (c{(fc), q(fe)), e'^(fe) = J^j^-'^'^H^ij) and A(fc) = ^ . e"*^'^^ A(j). We fur- 
ther require tjj = tjj. Then we have e^(fc) = e^{—k) = e(fc). Generically various realizations 
of order parameters are determined by the self-consistent equation. Here let us assume a 
form of the order parameters to discuss topological quantum phase transitions. 

The "spin" Hall conductance of the superconducting state on a lattice is given by the 
generalized TKNN formula as a = where C = ^ Jrj^2dSk ■ rot^ A^, = (fc|Vfcfc) 
P, 1^, D, |TT], |12|, |13[. The Chern number C, is given by a total vorticity of the quasiparticle 



Bloch function \k) for the negative energy quasiparticle, h{k)\k) = —E{k)\k), E{k) = 
Y^e(fcp~+~|A(fcyp since the quantum mechanical average is taken over the grand canonical 
ensemble. The integration is over the Brillouin zone which is topologically a torus T^, and 
dS = {dkxi 0, 0) X (0, dky^ 0) = (0, 0, dk^dky) is an infinitesimal area. The Chern number has 
an apparent topological expression which is written as a winding number of a relative phase 
of the Bloch function [0, || 

C = -N^inding = -y^^ <f dk- Vim log i^. 

V 271- JdR, \k)2 



I k) 

The summation is over the poles of (zeros of the the second component of the Bloch 
function 1^)2) in the Brillouin zone {Ri is a small area around each pole) ||^, P). 
By the Berry's parametrization[Q, the hamiltonian h{k) can be written as 

h{k) = R{k) ■ a 

where cr = (a^, (Jy, Cz) are the Pauli matrices and R = {R^, Ry, Rz) = 
(Re A(fc), — Im A(A;), e{k)). Then the expression for the Chern number is rewritten by these 
coordinates in R space as 

C = / dSn ■ TotR Ar 

where Ar = {R\'VrR) and \k) = \R{k)). In a particular gauge ( a phase convention of 

/ sin I \ 

the vector), it is written as |-R) = where {R,9,(f)) is a polar coordinate of R. 

\e'^ cos I / 

The integral region R{T'^) is a closed surface in three dimensional parameter space mapped 
from the Brillouin zone T^. As is well known, the corresponding vector potential defines a 
magnetic monopole at the origin in the three dimensional space as divrot Ar = — 27ri 6r{R) 



21] , I24 . Therefore, by the Gauss' theorem, we have another expression for the Chern 



number as 

C = [ dVSR{R) = -NiR{T^), O) = -N,,,,ring 
Ji?(T2) 

where Ncovering = N{R(T'^),0) is a degree of covering by the closed surface R(T^) around 
the origin O. This is another topological expression characterizing the negative energy band. 

Now let us calculate the winding number N^mding by this gauge. In the present gauge, 
Vfclm log = -Vfc0 and we have N^inding = ^Y.p_ §R(dR,) ' ^^^^^ 
of the second component of the Bloch function are defined by cosd/2 = 0, we rewrite the 
above equation as 



C* -^winding ^ ^ I (^R(T 



where is a negative 2;-axes in the R space (the Dirac string) and P_'s are intersection 
points of the closed surface R{T'^) and the Dirac string Z-. The integer I{R{T'^), z^) is an 
intersection number which gives the multiplicity of the local covering at P_'s. It is a positive 
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integer when the Dirac string is parallel to the normal vector of the surface x |^ 
and a negative integer when it is anti-parallel. N binding is apparently integral and zero if the 
origin is outside the surface. It is one of the advantages of the present expression. 
Now we have established a relation for the two topological integers as 

^covering ^winding 

N{R{T'),0) = Y,I{R{T'),z.). 

P- 

This relation is physically apparent but it confirms the consistency of the discussion. 

Demonstration As is discussed above, the Chern number is a topological quantity and is 
stable against small perturbation. Each superconducting state is characterized by an integer 
and the topological quantum transition is associated with topological change of geometrical 
objects, one of which is the relative phases of the Bloch functions of the quasiparticles and 
the others are the closed surface i?(T^) and the Dirac string z^. The topological transition 
of the superconductivity is clearly observed by investigating these topological objects. 

Now let us take following two examples of the order parameters to demonstrate the 
topological transition. 

Case I : singlet, Ajj = A^j, 

^i,i+x+y ^i+x,i-\-y i^xy 

Case II : triplet, Ajj — — Ajj, 

Ajji+x — '^^i,i+y — Aj; 

where A^2_y2 and A^y are real. The hopping f^j — t, is also real and only nonzero on the 
nearest neighbor links. We have e(fc) = —2t{coskx + cosky) — ji, for the both cases, and 
A(fc) = Ao + 2Ax'i-y'i{coskx — cosky) + 2iAxy{cos{k3; + ky) — cos{kx — ky)) for the singlet 
case (Case I) and A(fe) = 2iAa;(sinA;^ -\-ism.ky) for the triplet case (Case II). 

In Fig. 1, the relative phase of the Bloch function, Arg for the singlet case is shown. 

\K)2 

For small Aq, one can see the charge two vortices at the pole of |^ (circle) spht into 

\K2) 

two independent ones and contribute to the non trivial Chern numbers. Note that only the 
poles ( circles) of |^ contribute to the Chern numbers although the zeros (squares) also give 

\K)2 

vortices. When Aq is sufficiently large, they finally annihilate in pairs with negative vortices 



defined by the zeros of -rrr (squares). Anotlier topological objects, the closed surface R{T'^) 
and the Dirac string is shown for the singlet case in Fig.2. One can see the covering 
degree of mapping — > i?(T^) is two if Aq is sufficiently small. The intersection number 
of the closed surface i?(T^) and the Dirac string Z- is also two which can be easily checked 
for this example. In general, the order parameter may have more complicated momentum 
dependence, and calculating the covering degree of mapping is not trivial at all. Therefore 
the intersection number I{R{T'^),z_) helps us to determine the Chern number concretely. 
It is also easy to see that the Chern number vanishes if |Ao| or |yu| is sufficiently large since 
the surface i?(T^) does not include the origin. In this singlet case, the origin O, is doubly 
covered since R{k) = R{—k) ( A(fc) = A(— fc) ). It implies the selection rule AC = ±2[1T2|. 
Another interesting case is given by the A^y = 0. In this case, the energy gap collapses and 
the gapless Dirac dispersion prevent from determining the Chern number without ambiguity 
(if \fi\ is small). This situation is evident by the present discussion since the surface R{T'^) 
is collapsed into a diamond shaped two dimensional region R = XiVi + \2V2 — (0,0, /i), 
Ai,2 £ ["I)!]) "^1 = {Ax2_y2,0, —2t) and V2 = {—Ax2^y2,0,2t). When the origin O is on 
this region, the Chern number is ill-defined and the dispersion is gapless. This condition is 
clearly given by < 2t. Otherwise, there is an energy gap and the Chern number is zero. 

The triplet case, the case II, gives a little tricky situation. As shown in Fig. 3, the 
covering degree of mapping around the origin is +1 or —1 depending on the parameter detail. 
The surface R(T^) is self-intersecting and a local coordinate {e.j.{R) , ey{R) , n{R)) on the 
surface has a different chirality depending on the position of the surface where ea{R) = 
{a = x,y) are tangent vectors and n{R) is a normal vector of the surface at R. In other 
words, inside and outside of the surface are exchanged on intersecting lines. Of course, the 
surface i?(T^) itself is orientable always. At the transition, the selection rule for the triplet 
case is AC = ±1 generically. However, the transition in which the origin O passes through 
this intersecting point, say, changing the chemical potential n, the Chern number can change 
its sign as C : ±1 ^ ^1, (AC = ±2). 

In conclusion, we have demonstrated the topological quantum phase transition in super- 
conductivity on lattices by using several topological expressions for the Chern number. They 
are the winding number of the relative phase of the quasiparticle Bloch function N^inding 
and the covering degree Ncovering of mapping RiT"^) around the Dirac monopole at 

origin. The winding number also has another expression given by the intersection number 
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FIG. 1: Relative phases, Arg fr^ = Im log jyt' of the quasiparticle Bloch function in the Brillouin 

\K)2 \K)2 

zone ( solid square ). The circles and the squares include the poles and the zeros of of 
respectively. Note that only the circles ( poles ) contribute to the total Chern number. The spin 
singlet case with (a) Aq = 0, A^2_j^2 = t, A^y = 0.5t, /j, = (the vorticity of the pole at the origin 
is +2), (b) Aq = 2t, A^2_j^2 = t, Axy = 0.5t, /x = ( there are two vorticity +1 vortices at the 
poles). 

li^RiT"^), z^) of the closed surface R{T^) and the Dirac string Z-. 
We thank Y. Morita and Y. Seki for discussions. 
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FIG. 2: Mapped Brillouin zone R{T'^) by the Berry's parametrization. The monopole at the origin 
O and the Dirac string z- are also shown. The spin singlet case with Aq = 0, /S.^2_y2 = t, A^y = 
t, fi = —t. (a) is the total surface RiT"^), S [0, 27r], ky € [0, 27r] and (b)-(f) arc parts of the 
surface to show how it covers the origin. They are drawn for kx G [0, 27r] and (b): ky G [0, (c): 
ky e [f ,f ], (d): ky G [f ,f ], (e): ky G [f ,f ], (f): ky G [f ,27r]. The monopole O is doubly 
covered by surface R{T^). 

FIG. 3: Same as Fig.2 for the spin triplet case with Aq = 0, A-r = —t, jjl = —t. (a) is the total 
surface R{T'^), k^ G [0, 27r], ky G [0, 27r]. (b) and (c) are drawn for kx G [0, 27r] and (b): ky G [0, vr], 
(c): ky G [7r,27r]. 
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